AC-induced superfluidity 
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We argue that a system of ultracold bosonic atoms in a tilted optical lattice can become superfluid 
in response to resonant AC forcing. Among others, this allows one to prepare a Bose-Einstein 
condensate in a state associated with a negative effective mass. Our reasoning is backed by both 
exact numerical simulations for systems consisting of few particles, and by a theoretical approach 
based on Floquet-Fock states. 
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I. INTRODUCTION 



The study of ultracold atoms in optical lattices by now 
has opened up a promising new area of research on the 
borderlines between atomic and molecular physics, quan- 
tum optics, and condensed-matter physics: Many-body 
systems of paradigmatic importance to condensed-matter 
theory are being experimentally realized, with the preci- 
sion and flexibility offered by quantum optical tools, by 
means of interacting ultracold atoms subjected to pe- 
riodic potentials generated by standing waves of laser 
light. A hallmark example along this line is provided 
by the successful implementation of the Bose-Hubbard 
Hamiltonian [1, 2], and the observation of the predicted 
quantum phase transition from a superfluid to a Mott 
insulator 0, 0, [Bj]. Extrapolating these developments, 
it has been suggested to employ fermionic atoms con- 
fined by optical lattices in order to achieve a laboratory 
realization of the (fermionic) Hubbard model, and thus 
to tackle questions concerning high-T c -superconductivity 
by performing measurements on that system [6], while 
an exact numerical treatment of the very model lies still 
beyond existing computational resources. 



In this letter, we propose to investigate time-dependent 
many-body dynamics of ultracold bosonic atoms in a one- 
dimensional optical lattice subjected to both a static tilt 
and an additional time-periodic (AC) drive. We argue 
that such a system exhibits an effect which has no im- 
mediate solid-state analogue, namely, the appearance of 
a superfluid state in response to resonant forcing. We 
first discuss the underlying mechanism in a general man- 
ner, before presenting numerical calculations for small 
systems, and a quantum Floquet theoretical approach 
which both support our deductions. Finally we provide 
parameter estimates concerning the experimental verifi- 
cation of this effect. 



II. THE BASIC MECHANISM 



Ultracold bosonic atoms in a sufficiently deep one- 
dimensional optical lattice with M sites are described, to 
good approximation, by the Bose-Hubbard model p], Q 
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where b\ (b £ ) denotes the creation (annihilation) operator 
for an atom in the Wannier state located at the £th lattice 
site, hi = bgb £ is the number operator for that site, the 
hopping matrix element J > parametrises the strength 
of tunneling between adjacent sites, and U > quan- 
tifies the repulsion energy of a pair of atoms occupying 
the same site. For vanishing interaction strength, that is, 
for U / J —> 0, all N atoms condense into the lowest Bloch 
state of the lattice for temperature T — > 0, giving rise to a 

superfluid ground state |SF) = Y^ =1 b\j |0), 

where |0) denotes the vacuum state. On the other hand, 
for vanishing tunneling contact, U/J — > oo, the ground 

state is the Mott insulating state |MI) = ]\f =1 ^=f|0), 
assuming integer filling n = N/M . In the limit of a chain 
consisting of infinitely many sites, M — > oo, a sharp tran- 
sition between the superfluid regime with (quasi) long- 
range order and the Mott-insulating regime occurs at a 
critical value (U / J) c , accompanied by the emergence of 
a finite gap between the energy of the ground state and 
those of the excited states For a filling factor of one 
particle per site, N/M = 1, one finds (U/J) c ^ 3.4 Q. 
Remarkably, even calculations performed for quite small 
systems show a precursor of this quantum phase tran- 
sition: In Fig. [1] we display the energy eigenvalues for a 
system with N = M = 5, as functions of the scaled inter- 
action strength U/J. There is, of course, no sharp transi- 
tion in such a small system, but the gradual splitting-off 
of the ground-state energy is clearly discernible. 

When the Bose-Hubbard system (pQ) with M — > oo 
is subjected to both a static tilt amounting to an en- 
ergy difference per site of Kq, and to an AC force of 
frequency oo and amplitude K Ul the total Hamiltonian 
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FIG. 1: Energy spectrum versus interaction strength U/J 
for a small, untilted and undriven Bose-Hubbard-system (pQ) 
with N = 5 particles on M = 5 sites. The splitting-off of the 
ground state around U / J « 4 is a precursor of the quantum 
phase transition that occurs in an infinitely large system. 



becomes H(t) = Hq + H\{t), where 



A I 



ffi(t) = (i^o + i^u, cos(cjt)) J^ra* 



(2) 



Assuming first = 0, the static force effectuates a split- 
ting of the single-particle Bloch band into a sequence of 
single-particle states with equidistant energies E% = K§t 
forming the so-called Wannier-Stark ladder [8]. The 
corresponding single-particle orbitals are associated with 
creation operators c\ = ^ i ^i-i{2J / Ko)b\, where J a (z) 
is a Bessel function of order a. Since \Jq(2J/ Kq) 



l-2(J/K ) 2 + O[(J/K ) 4 } for J/K <C 1, one has cj w 6 
when i^o ^> «/, so that the Wannier-Stark states coincide 
approximately with the original Wannier states, indicat- 
ing that a strong tilt effectively destroys the tunneling 
contact between adjacent lattice sites. 

The principle we are exploiting in the present proposal 
relies on the fact that a tunneling contact disabled by a 
uniform tilt can be partially restored when the system is 
driven resonantly, i.e., when the driving frequency uo in 
eq. ([2]) is chosen such that 



Kq ~ vfajJ 



(3) 



with integer v, so that the energy of v quanta ftw 
bridges the energy difference between adjacent rungs of 
the Wannier-Stark ladder [9]. In that case the driven 
Bose-Hubbard system behaves approximately, in a time- 
averaged sense, like an undriven system with a modified 
hopping matrix element 



eff 



(-1)^(^/^)7, 



(4) 



provided the frequency is sufficiently high, so that both 
frjjj ^> U und frjjj ^> J. While such a modification of the 
tunneling contact is reminiscent of the AC-Josephson- 
effect [To|, 0j] , here the many-body interaction has a re- 
markable consequence: When \U/J G s\ becomes smaller 
than (U/J) c , the system becomes superfluid despite the 
presence of the strong tilt. Observe that this scenario 



also includes pure AC-forcing with Kq = as a special 
case for v = 0. For this case, which has been considered 
theoretically before [H|,[l3j], the expected Jo -modification 
of the hopping matrix element has been beautifully con- 
firmed in a recent experiment [14|, thus demonstrating 
the feasibility of exposing a Bose-Einstein condensate in 
an optical lattice to strong AC forcing without destroying 
its phase coherence. 

In order to verify the existence of AC-induced super- 
fluidity, we propose the following experimental protocol: 
(i) Initially the ultracold atoms are prepared in the su- 
perfluid ground state of a shallow one-dimensional op- 
tical lattice. Then the interaction parameter U/J is 
ramped up by increasing the depth of the lattice. At the 
critical parameter (U/J) c the system enters the Mott- 
insulator regime, and the particles localise at the lattice 
sites. For large values of U/J, the system approximately 
reaches the extreme Mott state |MI), again assuming in- 
teger filling n = N/M. (ii) Next, a strong static tilt 
is applied, e.g. by accelerating the lattice. Then U/J is 
switched down again to a small value, such that the sys- 
tem would fall into the superfluid regime if there were no 
tilt. However, it actually remains in a state close to |MI), 
which, in its turn, still remains an approximate station- 
ary state of the Hamiltonian, due to the strong localisa- 
tion of the Wannier-Stark states. (Hi) Finally the reso- 
nant AC force is turned on. In order to guide the system 
adiabatically into the AC-induced superfluid state cor- 
responding to the ground state of the effective undriven 
Bose-Hubbard Hamiltonian with the modified hopping 
element (j4j), that force should be turned on smoothly, 
starting from K u /Huj = 0. 



III. EXACT SIMULATION OF SMALL 
SYSTEMS 

We have computed the exact time evolution of a small 
system of TV = 9 particles on M = 9 lattice sites, yield- 
ing the filling factor N/M = 1. To analyze the time- 
dependent wave function \ip(t)), we calculate the single- 
particle quasimomentum distribution 



= M^ 6XP 



(i-j)p' 



h/c 



{m\i\b 3 \m) 



(5) 



where a is the lattice constant, and record it versus 
pa/(27rh) at times t that are integer multiples of the driv- 
ing period T = 2tt/uj. This distribution is 27r/z/a-periodic 
with respect to p; together with an envelope provided by 
the momentum distribution of the corresponding Wan- 
nier state it yields the momentum distribution of the 
system that can be measured experimentally by time-of- 
flight absorption imaging [1, 0, [5[ . It allows one to clearly 
distinguish between the superfluid and the Mott state: 
Whereas the momentum distribution of a Mott state is 
rather structureless, a superfluid state with (quasi) long- 
range phase coherence is characterised by sharp peaks 
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FIG. 2: Upper row: Three Brillouin zones of the quasimomentum distribution p(p,t), recorded at integer t/T, while the driving 
amplitude K u /huj increases linearly from to 10 between t = 0T and t = 200 T. Parameters for t > 0T are U/J — 1 and 
Huj/J = 15; N = M = 9. For z/ = 0, the system starts in a superfluid state, corresponding to a sharply peaked distribution. 
When Kuj/Huj passes a zero of Jo, the maxima of the distribution switch from the centers to the edges of the Brillouin zones, or 
vice versa. For v = 1 and v — 2, the system initially is prepared in a Mott-like state. Nonetheless, the sharp peaks emerging at 
later times, when | J e ff| is sufficiently large, signal the appearance of AC-induced superfluidity. Again, the peaks switch between 
the zone centers and the edges at the zeros of the Bessel function Ji or J2, respectively. — Lower row: Exact instantaneous 
quasienergies for the same situations, computed with N — M — 5. Observe that the changes in the momentum distributions 
above are precisely reflected by the band collapses. 



at the minima of the effective single-particle dispersion 
relation E e &(p) = —2J e ^cos(pa/h). 

Fixing the time scale such that the AC force with high 
frequency hco/J = 15 is turned on at t = 0, the system 
is initialised in the ground state of an undriven, untilted 
lattice with large interaction parameter U j J = 50 at time 
t = — 20 T. Simultaneously the lattice is tilted abruptly 
by K = vtvu, with v = 1 or 2. At t = -10 T the inter- 
action strength is decreased abruptly to U / J = 1. These 
preparatory steps lead to no significant change of the ini- 
tial Mott state. That state coincides to good accuracy 
(overlap « 0.95 or larger) with the ground state |MI) of 
the effective system for K^/hco = 0, i.e., for J e R — 0. 
For comparison, we also show results obtained for an un- 
tilted lattice (y = 0). In that case, the system simply 
is initialised at t = in its superfluid ground state for 
U/J=l. 

During the time intervall from 0T to 200 T the di- 
mensionless driving amplitude K u /huj = z is ramped up 
linearly from to 10. The upper row in Fig. [2] shows the 
resulting quasimomentum distributions ([5]), as obtained 
from numerical solutions of the time-dependent many- 
body Schrodinger equation. For v — 0, the system starts 
in a superfluid state. Accordingly, p(p, t) displayed in 
Fig. 2(a)| shows a high contrast right from the outset at 
t = 0, with sharp maxima at the Brillouin zone center, 



pa/(2,7rh) = 0, ±1, ±2, .... However, when the gradually 
rising amplitude K u (t)/huj passes the first zero of Jo (2) 
at z w 2.4, there is a short interval during which |J e ff| 
is so small, or \U / J e ^\ so large, that the system becomes 
Mott-like [12[- When the superfluid-like state reemerges 
at later t, the maxima of p(p, t) have shifted to the edges 
of the Brillouin zone. This finding is fully in accordance 
with the fact that J e ff is negative between the first and 
second zero of Jo, so that the minima of E e s(p) occur at 
the zone edges. When the driving amplitude rises even 
beyond the second zero of Jo (2) at z w 5.5, the original 
conditions are restored. We remark that this change of 
sign of the effective hopping matrix element, which has 
been clearly observed in the experiment [14], is a nontriv- 
ial consequence of the "dressing" of the system achieved 
through the time-periodic forcing: In effect, this amounts 
to placing the condensate in a state with a negative ef- 
fective mass, the sign of that mass being controlled by 
the forcing strength. 

When v = 1, the system initially is Mott-like at 
t = 0T. However, as soon as K^ii^/hw becomes suf- 
ficently large, the tunneling contact is restored, and the 
system beomes super fluid-l ike, as wittnessed by the sharp 
feature visible in Fig. |2(b)| This is precisely the effect we 
are focussing on. Observe that here J e ff is negative be- 
tween z = and the first zero of Ji(z) at z « 3.8, so 
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FIG. 3: Momentum distribution for parameters and protocol 
as in Fig. 2(c) but for an initial state with half filling, N — 



6 particles on M — 12 sites. Shortly after t — 100 T, when J e ff 
changes its sign, the signature of superfluidity is lost, since the 
system is unable to follow the rapidly changing ground state. 



that the negative- mass picture now holds even for small 
driving amplitudes. The results for v = 2 in Fig. |2(c) 
likewise follow the expected pattern. Thus, the com- 
putations summarised in the upper row of Fig. [2] can 
be explained in terms of an effective time-independent 
Bose-Hubbard Hamiltonian (pQ) equipped with the mod- 
ified hopping matrix element (|4|) at any instantaneous 
value of the driving amplitude, in combination with ap- 
proximate adiabatic following of the system's wave func- 
tion in response to the actual slow rise of that amplitude. 
This intuitive picture has a striking implication: A Bose- 
Hubbard system at half filling remains superfluid even 
for J/U ^ [1]. But when the sign of J e ff changes, the 
superfluid ground state changes abruptly, rendering adi- 
abatic following impossible. This is illustrated in Fig. [3l 
where we have plotted the momentum distribution for 
the same p arame ters and protocol as previously consid- 
ered in Fig. 2(c)| but for N — 6 and M = 12: For times 
slightly larger than 100 T, where J e ff becomes negative, 
superfluidity is lost. This expected outcome signals that 
the high degree of adiabaticity achieved in Fig. [2] for unit 
filling actually is a nontrivial many-body effect; it hinges 
on the presence of a Mott-like state when J e f?/U is close 
to zero. 



IV. QUANTUM FLOQUET THEORETICAL 
APPROACH 



We now give a formal explanation of the approxi- 
mate AC-induced modification of the hopping matrix el- 
ement (|4]). We separate the full Hamiltonian Ho + Hi(t) 
into a part which is site-diagonal, 



Hint(*) 



U 



J2n e (n, -!)+#!(*) , 



(6) 



and a part which describes tunneling between neighbour- 
ing sites, 



M-l 

-J £ (b\b 

£=1 



(7) 



We do not require that the resonance condition ([3|) is met 
exactly, which would be close to impossible in a labora- 
tory experiment, but rather admit a slight detuning, and 
choose the integer v such that \Kq — uhuu\ is minimised. 
Then one has 



K = vhuo + AK 
with |Aifo| <^ Now the set of Fock states 



KM) = E 4=10) 



me. 



(8) 



(9) 



where {n^} runs over all admissible sets of occupation 
numbers of the lattice sites, forms a basis for the diago- 
nalisation of the time-independent Bose-Hubbard Hamil- 
tonian (pp). However, when the forcing ([2]) is turned 
on with constant amplitude K^, the full Hamiltonian 
H(t) = H[nt(t) + ^tim is T-periodic in time, with T = 
2tt/uj. In that case the system's Floquet states, i.e., the 
complete set of explicitly time-dependent, T-periodic so- 
lutions \u(t))) to the eigenvalue equation 



H(t) - m \u(t))) = e \u(t))) , 



(10) 



take over the role of the stationary states; the wave 
functions \^{t)) = \u(t))) exp(— iet/h) solve the time- 
dependent Schrodinger equation. The spectrum of their 
quasienergies e is obtained by diagonalising the operator 
H(t) — ihd t in an extended Hilbert space of T-periodic 
functions, endowed with the scalar product 



((•I-)) 



dt(.\.) 



(ii) 



which combines the usual scalar product ( • | • ) with time- 
avaraging [15]. Accordingly, we employ a set of Floquet- 
Fock states 




KM) 



£ri£ + imcot 



(12) 



which diagonalise the part Hmt(t) — ihdt of the 
quasienergy operator: By construction, one has 



(({ne}, m\ H int (t) - \hd t |{n^}, m)) 



- rnhu . 



(13) 
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Observe that the quasienergy spectrum is periodic in e 
with period ftw; the integer m = 0, ±1, ±2, . . . introduced 
as a Fourier index in the basis states ([T2]) serves to distin- 
guish different Brillouin zones of that spectrum. Employ- 
ing the identity e lzsinip = ^2t=-oo elkip ^ k ( z ) > one easn Y 
obtains 

{{{n[},m'\b\b i+1 \{n e },m)) 

= (We}\ tyh+i KM) J s ( m -m0-^(^/M ( 14 ) 

with s = X^( n ^ — n i) = +1? since PJ> £ , 1 transfers one 
particle by one site to the left. Note that when taking the 
adjoint matrix element, both m, m' and {n^}, {n f £ } inter- 
change their roles, so that the sign s changes to —1 and 
the numerical value of the element remains unchanged. 
Hence, with respect to the index m the matrix of the 
quasienergy operator has a simple block structure: 

((K},m'|ff(t)-ifift|{^},m» 

({r4}| iJ e ff + mtvu \ {nt}) 

+ (l-^,m)(K}|y|{n,}) • (15) 

The diagonal blocks (m = m') are just the matrices 
of a Bose-Hubbard system in the standard Fock ba- 
sis (|9|), without AC forcing, but with modified hopping 
strength Q and possibly a slight residual tilt, as de- 
scribed by the effective Hamiltonian 

M-l jj M 

£=1 £=1 

+ AK J2£hi. (16) 

£ 

These blocks are shifted against each other by multi- 
ples of the "photon" energy Huj, as corresponding to 
the Brillouin-zone structure of the quasienergy spectrum. 
They are coupled by off-diagonal blocks provided by an 
operator V which differs from the tunneling term ([7j) such 
that the hopping element J is multiplied by a Bessel func- 
tion Jm^m'-v^u/huj) for each particle transfer directed 
to the left, as described by a combination b\b i+1 , and by 
^m'^m-yiK^lhbj) for each transfer directed to the right, 
corresponding to b\ +1 b^. 

The viewpoint advocated before, namely, the descrip- 
tion of the forced system Ho + H\ (t) = H[ nt (t) + Htun in 
terms of the effective time-independent Hamiltonian ([T6]h 
now amounts to neglecting all off-diagonal coupling 
blocks with m! ^ m. This is a viable approximation 
in the high-frequency regime, were (besides |AXo|) both 
relevant energy scales J and U are small compared to 
hco. Intuitively speaking, this condition guarantees that 
the AC drive is off-resonant with respect to low-order 
transitions. As expressed by the scalar product JTTJ), the 
replacement of the original Hamiltonian by H e R essen- 
tially relies on the averaging principle. 

To substantiate our reasoning, the lower row of Fig. [2] 
shows exact quasienergy spectra, obained from numeri- 
cal solutions of the eigenvalue equation ([TO]) with N = 



M = 5, for the same cases as considered in the dynami- 
cal simulations displayed in the upper row of that figure. 
These spectra are well matched to the effective Hamil- 
tonian ([T6]) with AKo = 0; in particular, the observed 
band collapses are accurately located at the zeros of the 
respective Bessel function J u . In the vicinity of these 
points the quasienergy of the effective ground state sep- 
arates markedly from the bands of excited states, thus 
indicating the presence of a Mott-like regime. The width 
of the quasienergy bands reflects the strength of the ef- 
fective tunneling contact: For v = 0, there is no tilt and 
the contact provided by Ht UYl is active for small K u /Huj, 
resulting in wide bands. In contrast, for v — 1 and v = 2 
tunneling is hindered by the static tilt, so that super- 
fluidity is enabled by a process akin to photon-assisted 
tunneling only when the AC drive is sufficiently strong. 

Moreover, H e ^ implies a simple scaling property: 
Quasienergy spectra obtained for different v and dif- 
ferent K^/hx), but corresponding to the same value of 
| J e ff |, should be almost identical. This fact is confirmed 
in Fig. [H where we have plotted spectra for v — 0, 1, 
and 2, and driving amplitudes K^/hw — 1.811, 0.711 
and 1.915, respectively, chosen such that | J e ff| = J/3 in 
each case. Indeed, the three spectra look strikingly alike; 
as expected; their gross features are obtained from the 
energy spectrum shown in Fig. [T] by rescaling the axes 
of that figure by a factor of 1/3, and then taking the 
eigenvalues modulo Hlu to account for the quasienergy 
Brillouin zones. 

Besides the existence of many-body Floquet states, a 
second key ingredient to our proposal is adiabatic follow- 
ing. In fact, Floquet states respect an adiabatic prin- 
ciple [3~6( | : the initial energy eigenstate evolves into the 
"connected" Floquet state when the AC drive is turned 
on sufficiently slowly. However, in the present scenario 
adiabatic following is endangered by the couplings asso- 
ciated with V, which we have neglected when retaining 
only the diagonal blocks of the quasienergy matrix ([T5]) . 
These couplings lead to many weak resonances among the 
instantaneous quasienergy levels. Hence, what appears 
as effectively adiabatic motion actually should be viewed 
as highly intricate Landau- Zener dynamics at multiple 
avoided level crossings, so that one may even deteriorate 
the quality of adiabatic following by reducing the rate 
of the parameter change in some cases. In the examples 
shown in Fig. [21 the squared overlap \(ip e ff \ 1 P{^))\ 2 °f the 
true wave function \ip(t)) with the ground state |^eff) of 
the instantaneous H e ^ (with AKo = 0) at integer t/T 
always remains larger than 0.26 w 0.89^, 0.64 w 0.95^, 
and 0.30 w 0.87^ for v = 0, 1, and 2, respectively, and 
N = 9. The extent to which this adiabatic principle can 
be exploited to guide the evolution of ultracold atoms 
in AC-driven optical lattices into certain desired target 
states is a subject requiring further studies, both exper- 
imental and theoretical. 
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(a)i/ = and K u /Huj = 1.811 (b)i/ = 1 and = 0.711 (c)i/ = 2 and = 1.915 



FIG. 4: Exact quasienergy spectra versus interaction strength U/J for a small system (N = M = 5) for Huj/J = 15, and 
Ko = ishu, with z/ = 0, 1, 2. Amplitudes K^/hcd are chosen such that | J e ff| = J/3. Observe that these spectra closely resemble 
the energy eigenvalues in Fig. [1] when the axes in that figure are rescaled by 1/3 and the eigenvalues are taken modulo Hlj. 



V. EXPERIMENTAL FEASIBILITY AND 
CONCLUSION 

In a laboratory experiment with optical lattices, some 
restrictions must be met in order not to leave the scope 
of the single-band Bose-Hubbard model: The static tilt 
must remain weak enough to prevent Zener tunneling on 
the time scale of the experiment, and both the frequency 
and the amplitude of the AC force must remain suffi- 
ciently low in order not to excite transitions to higher 
bands. Denoting the energy gap between the lowest 
two Bloch bands by E g , this translates into the require- 
ments U, J < K ,Huj,K u < E g , while U/J < (U/J) c . 
Considering a one-dimensional lattice with a depth of 
four atomic recoil energies and strong transversal confine- 
ment, as in the Zurich experiment |4(, one finds E g « 23 J 
and U/J ~ 614 a/d, where a denotes the s-wave scat- 
tering length of the atomic species employed, and d is 



the lattice constant. Thus, favourable conditions with 
U/J<1 are obtained when a/d < 0.0016. This demand 
appears severe, but not impossible to satisfy. 

In closing, we remark that the modification (|4|) of the 
hopping matrix element is very similar to the modifica- 
tion of atomic ^-factors by oscillating magnetic fields [ItJ , 
both mathematically and conceptually [18]. Understand- 
ing that modificaton has been crucial for developing the 
"dressed atom" -picture of atomic physics. Hence, the 
pioneering experiment [14j, having confirmed eq. (|4]) for 
v — with a Bose-Einstein condensate, eventually might 
lead to an analogous picture of "dressed matter waves" . 
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